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Abstract 

We give a new algorithm for merging sorted lists of monomials. To- 
gether with a projection technique we obtain a new complexity bound for 
the BM-algorithm. 

1 Introduction 

Vanishing ideals of points are of interest in many fields of mathematics — they 
are used in coding theory, in interpolation problems and even in statistics. Re- 
cently, the vanishing ideal of a set of affinc points has been studied in molecular 
biology |9j. 

The BM-algorithm [3] was proposed as a tool to make computations over 
vanishing ideals of points. When doing complexity studies of the BM-algorithm, 
one has to deal with arithmetic operations over the ground field and monomial 
manipulations. The number of arithmetic operations is reported [H |TU] to be 
proportional to nm'^, where n denotes the number of variables and m the num- 
ber of points. The number of integer comparisons needed for the monomial 
manipulations is reported [U [10] to be proportional to n^w? . In the biologi- 
cal applications, the coefficients of the points takes values in a finite field Zp 
and one usually has m n. Accordingly, one has begun searching for algo- 
rithms which are optimized for these situations. In [6 , an algorithm which 
uses 0{nm? -f m^) operations (arithmetic and integer operations arc treated 
the same) is given, while in , the same authors sharpened this bound so that 
it reads 0{nm^ + pm'^ +pm^\og{pm)) operations, where again, arithmetic and 
integer operations are treated the same. 

In this paper, we first make a thorough study of the complexity of comparing 
monomials. We restrict our analysis to the admissible monomial orders on n 
indeterminates given by invertible matrices with Z-coefficients. These orders 
associate to each monomial an n- vector of integers with the property that com- 
paring two monomials is the same as comparing the n-vectors lexicographically. 
Although this is a restriction on the set of admissible monomial orders, we re- 
mark that earlier complexity studies [U [TJ llOj have been performed only on a 
much smaller set of of admissible orders, e.g. lex, deglex and degrevlex. 

To give bounds for the monomial manipulations, we study algorithms for 
comparing lexicographically sorted n-vectors and give a fast algorithm for merg- 
ing lexicographically sorted lists of such vectors. Summation of ordered poly- 
nomials is one example of a situation where merge algorithms are used and 
thus, our merge algorithm could be used to speed up the computation of S- 



polynomials during the computations of a Grobner basis with respect to an 
ideal given by generators. 

The rest of the paper is focused on the BM-algorithm. We notice that the 
upper bound for the number of arithmetic operations during the BM-algorithm 
as given in J4] can be sharpened to read 

0{nm^ + min(TO, njirC'). 

After introducing a projection technique we show that the upper bound for the 
time complexity of the monomial manipulation part can be lowered to 

0(min(m, n)rn? log(m)) 

using one of the monomial orders lex, deglex or degrevlex. The factor log(m) 
comes from the model we use — that adding two integers bounded by m has 
time complexity 0(log(rn)). It follows that the bottleneck of the BM-algorithm 
are the arithmetic operations. 

Our method has better complexity than both the standard BM-algorithm 
and the methods optimized for the situations where m <C n. It is also more 
general than the methods in [HllTj since we do not assume that k is a finite field. 

Finally, we use our methods to give new complexity results for the FGLM- 
algorithm [1] and for the algorithms concerning ideals defined by functionals 
given in [TU| . 

Throughout the rest of the paper, we let S = k\xi, . . . , x„] denote the poly- 
nomial ring in n variables over a field k. The notion x", where a = (ai, . . . , q;„) 
will be used as short for x"^ • ■ • xJJ;". 

2 Monomial manipulations 

Comparisons of monomials is a crucial part in the complexity analysis of compu- 
tational algebra. However, different complexity models leads to different com- 
plexity results and this is problematic when it comes to comparing different 
algorithms. For instance, in [HIT], comparison of two monomials in n variables 
with respect to a monomial order is assumed to have time complexity 0(n), 
while in [TlllOj. the comparison of two monomials is assumed to be 0{n) integer 
comparisons. If one assumes an integer comparison to be made in constant time, 
then these complexities agree. But this is an unrealistic assumption, since the 
time needed for comparing integers is dependent on the size of the integers. 

After fixing a sound computational model, we will show that for a common 
class of monomial orders, it is possible to compare two monomials of degree 
bounded by m in time proportional to n max(log(m/n), 1). We then show that it 
is possible to merge two lexicographically sorted lists a and b of n-tuples in some 
set E with s and t elements respectively using at most min(s, t) -f n comparisons 
in E plus time proportional to max(s, t) log(n). This is better than the expected 
max(s,i)n E-comparisons. Combining these result, we obtain a method which 
merges two lists of monomials, sorted with respect to lex, deglex or degrevlex, 
in time proportional to min(s,i)rimax(log(m/n), 1) -|-max(s,i) log(n-l-m). 

2.1 Complexity model 

In computer algebra, it is often implicitly assumed that address and index arith- 
metics can be performed in constant time. This means that reading a byte from 
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any position in the memory is done in constant time, and reading k bytes is an 
0(fc)-operation. Such a model has advantages over a Turing machine, since it is 
easier to work with and even more reahstic in the cases when the memory on a 
modern computer is enough to handle the input data. We will use this model. 

However, when it comes to integer arithmetic (summation and compar- 
isons), the papers [TJ [21 21 [7] implicitly assume this to be done in constant 
time. But with the model used above, the correct complexities should read 
0(log(a) -t-log(6)) = 0(log(max(a, 6))) for summation and 0(log(min(a, 6)) for 
comparison for the positive integers a and b. For simplicity, we will assume 
that multiplication of the integers is 0(log(a) log(5)), although there are better 
bounds [5]. 

Inspired by jlOj . we will split the complexity studies in two parts. We will 
give arithmetic complexity for the arithmetic operations and time complexity 
for the monomial manipulations. We remark that the time complexity of per- 
forming / arithmetic operations always is at least 0{f). We do not deal with 
growth of coefhcients in the arithmetic operations, but refer the reader to [4]. In 
[T], techniques are discussed when k = Q using the Chinese remainder theorem. 

2.2 Monomial orders 

An admissible monomial order ^ on xi , . . . , x„ is a total order on the monomials 
which respects multiplication and has the unit 1 as the minimal element. A 
complete classification of admissible monomial orderings was first given in '11] . 
We will perform complexity analyses for a subclass of these orders, namely 
those defined by n x n-matrices of rational numbers which we define below. 
This is a proper restriction, since given an admissible monomial order and a 
natural number, there is an n x n-matrix of rational numbers that agrees with 
the given monomial order on all monomials whose degrees are bounded by the 
given number 8J. In the rest of the paper we will assume that all monomial 
orders are admissible. 

Let A ~ {dij) be an element in GL„(Q) with the property that the first 
nonzero entry in each column is positive. Then we can induce an order ~<a on 
the monomials in S by 

iff Aa* < A(3\ 

where Aa* < A/?* is the lexicographic order on Z", that is, {vi,...,Vn) < 
{wi, . . . , Wn) if vi = wi, . . . , Vi-i = Wi^i and Vi < wi for some i. Notice that 
multiplying a row of A by a positive integer does not change the order induced 
by A, hence we may assume that A is an integer matrix. 

Given A and a monomial x", we call the vector Aa^ the associated order 
vector to and we denote it by ov^(a;"). The simple observation 

ovAixix") = A{ai, . . . ,ai-i,ai + 1,q;j+i, . . .,«„)* = ovA(a;") -I- (a^, . . . ,a„i), 

gives us a handy formula for computing the order vector recursively. 

An important subclass of the orders defined above consists of (1) the lexico- 
graphical order, (2) the degree lexicographical order and (3) the degree reverse 
lexicographic order. These orders, called standard in the rest of the paper, are 
the common most used ones in computer algebra and computer algebra systems 
have them predefined. 
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To clarify the notion of a matrix representing an order, we show below the 
matrix representations for the standard orders on xi,X2,Xs, all of which satisfy 

Xi )^ X2 y X3. 




Alex - 1 ,Adeglex =1 ,A 





? -^degrevlex 



If xi :^ • • • >- a;„ for a standard order -<, we have oviex{x") = a, ovdegiexix") = 
(J2t ai, • • • , ctn-i) and ovdegreviex = iJ2t • • • ) -"2)- It is casily seen 

how to compute ov{xiX") given ov(a;") for these orders. 

In general, if Xi^ y ■ ■ ■ >- Xi^ and -< is standard, we will assume that the 
sequence ii, . . . ,i„ is given a priori and that .t" means a;"^ • • - a;"". We do not 
make this assumption when -< is given by a matrix. As indicated in the intro- 
duction, we will see that the complexity analysis of the monomial comparisons 
is dependent on whether the order is standard or not. 

When is a monomial order on xi, . . . , 3:„ it will be useful to restrict -< to 
a subset of the variables. If Ess is such a subset, we write ^ess to denote the 
restriction of -< to Ess. 



2.3 Comparing vectors of integers 

Since we assume that comparing two monomials is the same as lexicographically 
comparing their associated order vectors, we will now focus on comparing vectors 
of integers. 

To be able to prove the next lemma, recall that the number of bits needed 
to represent an integer a is 



numbits(a) = 



2 if a = 

[log2(|a|)J + 2 otherwise 

where one bit is used to represent sign. 



Lemma 2.1. Let a = (ai, . . . , Q!„) be a vector of integers ai. Let m = J2i 
Then 

E, . . . ^ f 3n if m<2n 

^ ^1 — y n log2 (m/n) + 2n otherwise 

Proof. Suppose that a contains k non-zero entries. Without loss of generality, 
we may assume that ai=Q \ii> k and ^ if i < fc. We get 

k n 

numbits(aj) = numbits(aj) -|- numbits(O) 

j j=l j=k+l 

k k 

= I^LlogsKOJ + 2A; + 2(n - A;) < ^ logaK) + 2n. 
i=i j=i 

Now ^j-i log2(Q;j) = log2(ai • • • ccfe) and since 

ai • • • afe < m/k- ■ ■ m/k, 
•• V ' 

k times 
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we conclude that 

numbits(aj) < /(fc), 

3 

where 

/(fc) = fclog2(m/fc) + 2n. 

We see that f'{k) = for k = m/2 and that /(m/2) is a maxnTiuni. If m < 2n, 
then ^^(numbits(Q;i)) < /(m/2) = m/2 + 2n < 3n. If 2n < ni, then observe 
that f{k) is a monotone increasing function on the interval [l,m/2] so that 
^•(numbits(Q;i)) < f{n) — nlog2(m/n) + 2n. □ 

We can formulate Lemma 12.11 in a more compact way by saying that the 
number of bits needed is proportional to nmax(log(m/n), 1). 

Lemma 2.2. Let v = (wi, . . . , w„) and w = (wi, . . . , m„) be vectors of integers 
such that \vi\ — mi and \wi\ — m2. Let m = max(mi,m2). With time 
complexity 0(nmax(log(m/n), 1)) we can determine if v and w differ and if 
they do, we can determine the index i where they differ. 

Proof. We proceed as follows. Compare vi and wi. If they differ, we stop and 
return the index. Otherwise, continue until we reach an index i such that u,; ^ Wi 
OT V — w. Let us now analyze the complexity of this procedure. To compare 
the vectors w an w, we need to compare at most all entries, that is, we get time 
complexity proportional to ^j(numbits(i'i) + numbits(u'i)). We rewrite this 
sum as numbits(ui) + numbits(wi) and use Lemma [2.11 and the remark 
thereafter to conclude that this sum is dominated by n(max(log(mi/n), 1) + 
max(log(m2/n), 1)). 

□ 

2.4 Comparing monomials with respect to a standard or- 
der 

Lemma 2.3. Suppose that -< is a standard order. An upper bound for computing 
ov(xia;") is 0{n\og{rn)), where m is the degree of x°' . 

Proof. To compute ov(x") is 0(1) when ~< is lex, while for a degree order, we 
need to sum all entries in order to compute m. Since all entries are bounded by 
m, the lemma follows. 

□ 

Lemma 2.4. Suppose that -< is a standard order. An upper bound for computing 
ov(xiX°') given ov(a;") is 0(log(m)), where m is the degree of x" . 

Proof. We need to increment at most two entries (in case of a degree order). 
The lemma follows since all entries are bounded by m. □ 

Lemma 2.5. An upper bound for the time needed to compare two monomi- 
als a;" and given oy^{x°') and o\-^{x^) with respect to a standard order is 
0(nmax(log(m/n), 1)), where m — max(^j- a^, /3i). 

Proof. Immediately from Lemma [2.2l when -< is lex. When -< is one of the degree 
orders, the first entry equals the degree and the sum of the rest of the entries is 
bounded by m, thus we get the same complexity in this case as well. □ 
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Lemma 2.6. Let ^ be a standard monomial order and let Ess — {xj-^ , ■ • ■ , Xj—} 
be a subset of the variables such that Xj-^ ■ ■ ■ >- Xj—. To compute ~<ess is 0{\). 



Proof. Immediate. 



□ 



2.5 Comparing monomials with respect to a matrix order 

Lemma 2.7. An upper bound for the time needed for determining oya{x'^) 
by computing Aa* is proportional to max(log(TO/n), 1) log(c), where c — 
max(|aij|) and m ~ ai. 

Proof. To determine the i'th index of ov^(a;"), one needs to compute anai + 
• ■ ■+ai„a„. The cost for the multiplication is proportional to log2(c) log2(Q!i) 
— log2(c) log2(Q;i). By an argument similar to the proof of Lemma 12.11 
it follows that this expression is dominated by an expression proportional to 
rt max(log(r7i/n), 1) log(c). The cost for the addition is negligible and since i 
runs from 1 to n, the lemma follows. □ 

Lemma 2.8. An upper bound for the time needed for computing 0VA{xix") 
given ovAi^") is proportional to nlog{cm), where c — max(|a.y|) and m — 



Proof. From the recursion formulas given in section 12.21 we see that on a {x" ■ 
Xi) = ov^(a;") + (aii, . . . , a„i), that is, we need to do n summations of integers 



Lemma 2.9. An upper bound for the time needed to compare two monomials x" 
and x^ with respect to a matrix order defined by A given ov^(a;") and oya{x^) 
is proportional to nlog(cm), where c = max(|aij|) and m = ai. 



Let ^ be a monomial order given by a matrix and let Varord_;(n) be the 
cost of determining Ji, . . . , z„ such that Xi-^ )^ ■ ■ ■ >~ Xi^. When -< is a standard 
order we assume that zi, . . . , z„ was given as input, so that Varord^(n) is 0(1). 
However, we do not assume this for a general order given by a matrix. Instead 
we have 

Lemma 2.10. For a general ordering given by a matrix, we can compute 
Varord^(?i) in time 0{n'^log{c)\og{n)). 

Proof. To compare Xi -< Xj is the same as comparing the i'th and the j'th 
column of the matrix A defining An upper bound for the comparisons is thus 
0{nlog{nc/n) +n) = 0(nlog(c)). Since sorting is 0(nlog(n)) comparisons, the 
upper bound becomes 0{n^ log(c) log(7i)). □ 

Lemma 2.11. Let -< be a monomial order given by a matrix A and let Ess — 
{xjj^, . . . ^Xj^} be a subset of the variables such that Xj-^ )^ ■ ■ ■ >- Xj—. To deter- 
mine annx n-matrix Aess such that ^ess is given by A-^ss can be done using 
0{nrT^) arithmetic operations over Q. 



bounded by cm. 



□ 



Proof. We have n comparisons of integers bounded by cm. 



□ 
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Proof. Clearly the n x n matrix A obtained by keeping the columns ji, . . . , 
defines ^ess and it has rank n. Suppose that the j'th row of A can be written 
as a linear combination of the rows whose indices are less then i. Let a and 
/9 be two order vectors with respect to Ess. Then, if Aa and A(3 agree on the 
first z — 1 rows, then they also agree on the I'th row. Hence the z'th row is 
superfluous. Thus, to determine A^ss is the same as rowreducing A^ which has 
arithmetic complexity 0{rm^). □ 

2.6 Merging sorted lists of monomials 

Recall that the merging of two sorted lists a and h is done using 0(max(i, s)) 
comparisons, where t is the number of elements in a and s is the number of 
elements in 6. When a and 6 are lists of n-tuples of elements in some set and 
a and h are sorted lexicographically in increasing order, then we can improve 
the classical merge algorithm. We will use this result to analyze the cost of the 
monomial manipulations in the BM-algorithm. 

Let w = (ui, . . . , u„) and w — {wi, . . . , w„) be two n- vectors of non- negative 
integers, li v ^ w, let A(w, w) be the first index where v and w differ, li v = w 
let A('i;, w) = n + 1. 

Lemma 2.12. If u < v, u < w and A(u, w) < A(u, w), then u < v < w and 
A(t;, w) — A(u, w). 

Proof. Let k = A{u,w). Then ui = wi, . . . ,Uk-i = iffc-i and Uk < w^. Since 
k < A(w, v), we have ui = vi, . . . ,Uk — v^. Hence vi — wi, . . . , Vk-i = Wk-i 
and Vk < Wk and thus v < w and A(u, w) = k = A{u, w). □ 

Lemma 2.13. Ifu<v, u <w then A{v,w) > mm(A{u,v), A{u,w)). 

Proof. If A{u,w) < A{u,v) or A(u,v) < A(u, w), then the lemma follows by 
Lemma 12.121 Otherwise, A{u,v) = A{u,w) implies that v and w agree on the 
first A{u,v) positions. □ 

Lemma 2.14. Let a — (oi, . . . , at) be a list of n-tuples of elements in an ordered 
set Yj. Suppose that a is sorted lexicographically in increasing order and that we 
are given A{ai, ai^i) for i = 1, . . . ,t — 1. Let b be any element in E". Using 
0(t + n) comparisons of elements in S plus time proportional to ilog(n), we 
may find an index i, < i < t, such that ai < ■ ■ ■ < ai < b < a^+i < ■ ■ ■ < at 
and A(ai, b)when i > 1 and A(&, Oi+i) when i < t. 

When S is the set of non-negative integers and Vi < m for all v ^ a and 
^ ™7 upper bound for the time complexity is 0(nmax(log(m/n), 1) + 
t log(max(r7,, m))). 

Proof. The proof contains three parts. We (a) give an algorithm, (b) prove its 
correctness and (c) show that the complexity of the algorithm agrees with what 
was stated in the lemma. 

The algorithm 

At stage 0, compute k = A(ai, b). If ai < b, continue with stage 1. If 6 < ai we 
stop and return and A(&, ai) = k. 

At stage i, 1 < i < i, we suppose that Oi < b and that A(ai, b) is computed. 
If Oi = tti^i (which is equivalent to A(ai, Oi+i) = s + 1), we have a.; < a^+i < b. 
Thus, we set A(ai+i, b) = A{ai, b) and continue with stage i + 1. 
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Else, if A(ai,5) > A{ai, Ui+i) then ai < b < a^+i and A(6, Oj+i) = A{ai,ai+i) 
by Lemma 12.121 Thus, we stop and return i, A{ai,b) and A(6, a^+i) (which 
equals A(ai, a^+i)). 

Else, if A(ai,6) < A(aj,ai+i) then < a^+i < 6 and A(a,i+i,6) = A(ai,6), 
again by Lemma [2.121 We set A(aj+i,6) = A(aj,6) and continue with stage 
i + l. 

Else, if A(ai,6) = A(ai,ai+i) fc, compare bk and a^+i^fc, 6^+1 and a^+Lk+i 
and so on until we either conclude (1) bk> < ai^i^k' (2) bk' > fli+i.fe/ or (3) 
b = fli+i, where in cases (1) and (2), k' = A(6, a^+i). In case (1), we have 
b < ai+i. We stop and return i, A{ai,b) and A(ai+i,6) = fc'. In case (2), we 
have b > a^+i. We set A(ai+i,6) = k' and continue with stage i + l. In case 
(3), we have b — ai^i. We stop and return i, A(ai, b) and A(ai+i, 6) = t + 1. 

At stage t we have at < b. Thus, we stop and return t and A(at, b). 

The correctness of the algorithm 

By construction. 

Complexity of the algorithm 

There are two key indices that we update during the algorithm. The first (z) 
refers to a position in the list a, the second (fc) refers to a position in the vector 
b. Notice that after each comparison, either i or k is increasing. Since both i 
and k are non-decreasing it follows that the number of comparisons of elements 
in E is at most n + t. 

The number of A-comparisons during the algorithm is one per stage, that is, 
at most t. Every such comparison consists of comparing integers bounded by n. 
We conclude that the time needed for the integer comparisons is proportional 
to t\og{n). 

Let now S be the set of non-negative integers. Everytime we increase i, we 
make a comparison of an integer bounded by m, this gives time proportional 
to Hog(m). However, when increasing fc, we are in a situation where aik — 
bk- Hence, the total timed used for the increasings of k is proportional to 
n max(log(m/n), 1) bv l2.2l Only once during the algorithm we will compare aa- 
and bk to conclude that they differ, this cost is 0(log(7Ti)). It follows that an 
upper bound for the algorithm is proportional to 

t\og{n) -l-ilog(TO) -f- nmax(log(m/n), 1) + log(m) 

= t log(max(m, n)) + n max(log(m/n), 1). 

□ 

Theorem 2.15. Let a = (ai, . . . , at) and b — (bi, . . . , b^) be two lists of n-tuples 
of elements in an ordered set E. Suppose that a and b are sorted lexicographically 
with respect to the order < in E. Suppose that we are given A{ai, Oi+i) for 
i — 1, ... ,i — 1 and A(6i,5i+i) for i — 1, ... ,s — 1. We can merge a and b 
into a new list c and compute the sequence A(ci, C2), A(c2, C3), . . . using 0{sn + 
t) comparisons plus time complexity 0{tlog(n)). When E is the non-negative 
integers and Vi < m for allv G aUb, an upper bound for the time complexity 
of the algorithm is 0(min(s, i)n max(log(m/n), 1) + max(s, t) log(n + to)). 

Proof. Suppose, without loss of generality, that s < t. Let ii be the index 
returned after calling the algorithm in Lemma 12.141 with a and bi . Without 
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affecting the complexity of the algorithm, it is clear that we can modify it to 
give a list c = (oi, . . . , a^j^, 61). Suppose that i2 is the index returned after 
calling the algorithm in Lemma l2. 141 with Oij+i, . . . , a* and 62. Again, it is clear 
that it is possible to update the list c without affecting the complexity so that 
it reads (oi, . . . , a^^ , 5i, ajj^+i, . . . ,0^2,62). If we proceed in this way we obtain 
the sequence 

ai < • • • < flji < ^1 < a-i^+i < ■■ ■ <ai^ < b2 < ■ ■ ■ < bg < a^^+i < • • • < a*. 

Since the algorithm in Lemma |2. 141 returns A(6j, ai^.+i) we only need to check 
the case when a^^- = 0^^.,^^ in order to conclude that A(ci,C2), A(c2,C3), ... is 
computed as a side affect of the calls to the modified algorithm. But when 
o-ij — we have that bj and bj^i are consecutive and hence A(6j,6j+i) is 

already computed by assumption. 

Although it does not affect the complexity, stage of the algorithm given 
in Lemma 12.141 can be modified. When calling with a^^+i, . . . , and bj^i 
we can use that A(6j, ai^.+i) already is computed. Indeed, A(6j+i, Oi^+i) > 
min(A(6j, ai^.+i), A(&j, tij+i)) by Lemma |2.13[ so we could call the algorithm in 
Lemma 12.141 with the extra parameter min(A(&j, ai^.-|_i), A(6j , 6^+1)) to speed 
up the computation of A(6j+i, a^^+i) in stage 0. 

Since we make s calls to the algorithm in Lemma l2.14l (or to be more precise, 
to the modified algorithm as defined above), we make ((ii + 1) + n) + ((i2 + 1 — 
(ii + l)+n) + - • ■ + {{is + l — {is-i + l)) + n) = {is + l) + sn < t+sn comparisons of 
elements in E. The time complexity for the integer comparison part is by Lemma 

[SUlproportional to (ii + 1) log(n) + («2 + l-(«i + l)) log(7i)H h(is + l-(«s-i + 

1)) log(n) = {is + 1) log(n) < t log(n). If E is the non-negative integers, then the 
time complexity of the algorithm becomes 0(sn max(log(m/7i), l)+t log(7i+m)) 
by Lemma O □ 

We now give two applications of the new merge algorithm. The idea of the 
first example is to make the algorithm clear to the reader, while the second 
example shows the strength of the algorithm. 

Example 1. Suppose that we want to merge the lists 

/ 2 2 32 2 22 3\ 

a — (X1X3X4, XiX^, X1X4, X1X2X5, XiX2X4^X5, Xi) 

and 

b— {xi,xixl,xlx2XiX5,xlx2xlx5) 

of monomials, sorted in increasing order with respect to lex and xi > X2 > 
X3 > X4 > X5. Since oviexix") = a, we will use the exponent vectors. Thus, in 
accordance with the notation above, 

ai = (1,0, 2, 2,0), 02 = (1,0, 3, 0,0), 03 = (2, 0, 0, 1, 0), 04 = (2, 1, 0, 0, 1), 
as = (2, 1,0, 2,1), = (3, 0,0, 0,0) 

and 

61 = (1,0,0, 0,0), 62 = (1,0, 2, 0,0), 63 = (2, 1,0, 1,1), 64= (2,1,0,2,1). 

We begin by comparing bi and ai. We see that A(fei,ai) — 3 and bi < ai. 
Since A(&i,62) = 3, we compare 623 o-nd 013. They are equal, so we check the 
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fourth index and conclude that 62 < ai. Now A{b2,b^) — 1; so we conclude that 
0-1 < ^3 only by checking the first index. Since A(ai,a2) = 3 > 1, we also have 
a2 < 63. But A(a2, 03) = 1, so we compare 03 and 63 from index 1 and conclude 
that aa < 63 anii A (03, 63) = 2. So far we have 

&i < 62 < ai < a2 < aa- 

We see that A(a3,a4) = 2, thus we compare 632 and 042. T/iey are equal, so 
are 633 and 643, &wt 634 > a44, hence 04 < 63 and A(&3,a4) = 4. S'mce also 
A(a4, a^) = 4, we compare the fourth index ofb^ and 05 to conclude that b^ < a^. 
Finally, since A(63, 64) = 4, we check the fourth and fifth indices of 64 and 05 
to conclude that 64 = 05 . We have 

bi < b2 < ai < 02 < 03 < 04 < bj < b^ < 05 < og. 

and the sequence of differences is 

3,4,2,1,2,4,4,6,1. 

Example 2. Let n = 2s. Let f — xix^ + X2Xs and let g = X1X2 + X2X3 + 
■■■ + Xn-iXn. Let -< be degrevlex with respect to Xi >~ ■■■ >- x„. We see 
that the terms of f and g are written with respect to this order. Suppose that, 
during a Grobner basis computation, we want to compute the S-polynomial of 
f and g, that is, we want a sorted expression of S{f,g) ~ X2f — x^g. This is 
the same as merging x\xs and (— X2X§, — x|x4, —x^x^x^, —x^x^xq, —x^XqX-j, . . ., 
— a;3a;„_2a^n, — 2;3a;„_ia;„) (together with an arithmetic operation in the case of 
equality). For simplicity, we write these expressions as lists of order vectors and 
omit the coefficients. We get 

a= (3,0^_^,-l,0,...,0,™2) 

s times 

and 

6= ((3,0,..., 0,-2,-1), (3,0,..., 0,-1, -2,0), (3,0,..., 0,-1, -1,-1,0), 

(3,0,..., 0,-1, -1,0, -1,0),..., (3, -1,-1,0,..., 0,-1,0)). 

We assume that we are given the sequence of differences for g a priori. Since 
the sequence of differences is closed under multiplication with a monomial, we 
obtain the sequence of differences for b. It reads (n — 2, n — 3, . . . , 2). Using the 
algorithm in Lemma \2.14\ we first compare a and bi. After s + 2 comparisons, 
we see that a ^ bi. Since A{bi, 62) — n — 2 and A(a, 61) = s + 2, we get (after 
comparing n — 2 and s + 2) that a >- b2 and that A(a, 62) = s + 2. Continuing 
this way we see that a ^ 61, a &2; a ~< b^, . . . ,a ^ &s-3, and A(a, bi) = s + 2 
for i — 1, . . . , s — using 

s + 2 + 1 + ■ ■ ■ + 1 = n - 1 

s — 3 times 

comparisons. We have A(&s_3, 65-2) = s + 2 and since also A(a, bs-3) = s + 2, 
we compare 65-2,5+2 cmd as+2 to conclude that a >- 62-2- In total, we have used 
n— l + 2 = 2s+l comparisons. 

Proceeding in a naive way we would use (s + 2) • (s — 2) = — 4 comparisons. 
It should be remarked that one needs extra cost for the bookkeeping of A{a,bi), 
for i ~ 1, . . . , s — 2, which is of the same magnitude as the comparisons, that is, 
to be fair, we should compare 2s + 1 + s — 2 = 3s — 1 with s^ — 4. 
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3 The BM-algorithm revised 



We first fix some notation. If x°' is a monomial, we define its support, which 
we denote by supp(x"), to be the set of all Xi such that ai > 0. If M is a set 
of monomials, then we define the support of M to be the union of the supports 
of the elements in AI. Let / be an ideal in S and let B be any subset of S 
such that [B] ~ {[b] : & G B} is a vector space basis for S/I. Here [b] denotes 
the equivalence class in S/I containing b. If s is an element in S, its residue 
can be uniquely expressed as a linear combination of the elements in [B] , say 
[s] ~ "^Cilbi]. The S'-element ^Cibi is then called the normal form of s with 
respect to B and we write Nf (s, B) Cibi. We abuse notation and say that 
B (instead of [B]) is a basis for S/I. 

Let -< be an (admissible) monomial order. The initial ideal of /, denoted by 
in(/), is the monomial ideal consisting of all leading monomials of / with respect 
to -<. One of the characterizations of a set G being a Grobner basis of an ideal 
I with respect to a monomial order -< is that G C I and that the leading terms 
of G generate in(/). An old theorem by Macaulay states that the residues of 
the monomials outside in(/) form a k-basis for the quotient S/I. The set of 
monomials outside in(/) is closed under taking submonomials. A consequence 
of this is that we always have the unit 1 in such a basis (given dimiic(iS'//) > 0). 
Sets of monomials which are closed under taking submonomials is called order 
ideal of monomials, abbreviated by OIM. 

If p is a point in k" and / is an element of S, we denote by f{p) the evaluation 
of / at p. When P = {pi, . . . is a set of points, /(P) = (/(pi), . . ■ ,/(Pm))- 
If F = {/i, . . . , fs} is a set of elements in S, then F{P) is defined to be the 
s X m matrix whose i'th row is fi{P). 

The vanishing ideal I{P) is the ideal consisting of all elements in S which 
vanishes on all the points in P. If /i and /2 are two elements in S and [/i] — [/2] 
in S/I{P), then fi{p) — f2{p) for p £ P. That a set [B] of m elements is a 
k-basis for S/I is equivalent to dimiic(i3(P)) = m. 

The BM-algorithm takes as input a set of points in k" and a monomial 
order. It returns a Grobner basis G of I and the set B of monomials outside 
in(J). The BM-algorithm was first given in [3J. During the years it has been 
reformulated and modified. In the paper the ideas of the BM-algorithm 
was used to switch between different Grobner bases of a zero-dimensional ideal. 
In the unifying paper [TUj it was shown that both the BM- and the FGLM- 
algorithm can be seen as an algorithm that computes a Grobner basis from an 
ideal defined by functionals. The complexity studies given in [10' apply to the 
BM-algorithm and in fact, [lOj is by tradition the paper which one refers to 
when complexity issues of the BM-algorithm are discussed. However, most of 
the complexity studies in [lOj are done by referring to the paper f5J. 

We will first discuss the complexity studies of the BM-algorithm and post- 
pone the connection with ideals defined by functionals to section [3.31 

3.1 Two formulations of the BM-algorithm 

We first try to describe the ideas of the BM-algorithm without going into details. 
The algorithm uses a list L of possible basis-candidates, a list G of the partial 
Grobner basis and a list B of the partial basis-elements. In the main loop, one 
checks if 1{P) G span]]5(_B(P)), where / is the least element in L with respect to 
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If it was, write I as a linear combination "^^Cibi of the elements in B and 
insert I — Cibi into G. Otherwise, insert I into B and update L with new 
basis candidates. 

We give below the formulation of the BM-algorithm as given in When 
the algorithm terminates, G is the Grobner basis with respect to -< and B is 
the complement of the initial ideal with respect to -< . 

CI Start with empty lists G~B — R— []a list L — [1], and a matrix 
C = (cij) over k with m columns and initially zero number of rows. 

C2 li L = [], return the pair [G, B] and stop. Otherwise, choose the monomial 
t = min^(L), the smallest according to the ordering -<. Delete t from L. 

C3 Compute the evaluation vector {t{pi) , . . . , t{pm)) G k™, and reduce it 
against the rows of C to obtain 

(wi, . . . = . • ■,t{pm)) - ^ai(Qi, . . .,Cira) at e k. 

i 

C4 If (ui, . . . , Vm) = (0, . . . , 0), then append the polynomial t — a^ri to the 
list G, where is the i'th element of R. Continue with step C2. 1 

C5 If (wi, . . . , Vm) / (0, . . . , 0), then add (wi, . . . , Vm) as a new row to C, and 
t — aiTi as a new element to R. Append the power product t to B, and 
add to L those elements of {xit, . . . which are neither multiples of 

an element of L nor of in(G). Continue with step C2. 

The authors in [1] claims that this is the same as the algorithm restricted to 
the BM-situation which appeared in [10 , but this is not exactly the case. To 
get the algorithm given in [10^ restricted to the BM-situation using the five-step 
description, we need to reformulate steps 2 and 5. 

C2' If L = [ ], return the pair [G, S] and stop. Otherwise, choose the power 
product t = min_^(L), the smallest according to the ordering ^. Delete t 
from L. If t is a multiple of an element in in(G), then repeat this step. 
Else, continue with step C3. 

C5' If (wi, . . . , Vm) / (0, . . . , 0), then add (wi, . . . , Vm) as a new row to G, and 
t — aiVi as a new element to R. Append the power product t to B, and 
merge {xit, . . . , Xnt} and L. Continue with step C2'. 

Lemma 3.1. The output from the two algorithms given above agree and G will 
be a reduced Grobner basis for I(P) and [B] will be a basis for S/I{P). 

Proof. It was proved [T^ that using steps C2' and C5', G will be a reduced 
Grobner basis and [B] will be a basis for S/I. To prove that the two algorithms 
agree, one needs to show that the elements which are multiples of an element 
in L and not inserted into the list L during stage C5 are not needed in the 
computations. We now sketch this argument. Suppose that Xit is an element 
which is a multiple of an element in L during the algorithm using steps C2 
and C5 and suppose also that Xit is in the complement of in(/). Then all 
submonomials of Xit is in B. Suppose that u is the largest. Clearly |u| = |a;it| — 1. 
Thus, Xit = XjU for some Xj. At some stage during the algorithm using steps 
C2 and C5, u will be the least element in L and accordingly, no submonomials 
of Xit will be left in L. Thus, XjU = Xit will be added to L in step C5'. □ 
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These two variants give the same arithmetic complexity, which is reported 
[TnHHE] to be 0{nm^) arithmetic operations. We now give a better bound. 

Proposition 3.2. The arithmetic complexity of the algorithms given above agree 
and an upper bound is 0{nm? + min(m, n)m?). 

Proof. The arithmetic operations are performed in steps C3 and C4. Step C3 in- 
volves an evaluation and a row reduction. To compute . . . , t{pm)) requires 
m multiplications, since t — Xit' for some Xi and some t' and t'{pi), . . . ,t'{pm) 
already has been evaluated. The row reduction requires 0{m?) arithmetic op- 
erations. Step C4 consist of expressing a^ri in the basis and is also an 
0{m?) operation. Notice that the number of calls to C4 is exactly \B\ ~ m, 
the number of calls to C5/C5' is exactly \G\ and thus, the number of calls to 
C3 is exactly m -I- \G\. Thus, the number of arithmetic operations is propor- 
tional to (to + \G\) ■ m?), so it is enough to prove that \G\ is proportional to 
n + min(TO, n)m. 

The number of variables in B is at most min(r7,, m — 1) since the basis consist 
of TO elements and 1 G _B. Accordingly, if we let s denote the number of variables 
in in(G'), then n — min(n,TO — 1) < s <n {s = n only if to = 1). Now notice that 
if Xit is an element inserted into L, then, since Xi < Xit and Xi is inserted into L 
during the first step of the algorithm, Xi will be treated before Xit. This shows 
that when we add a new element t to B and {xit, . . . , Xnt} to L, we know that 
at most ri — s of these elements would be added to in(G), because if Xi G in(G'), 
then Xit G in(G'). 

Since we add to elements to B, we see that the elements of degree more than 
one in in(G) is at most m{n — s). We conclude that the number of elements in 
in(G) is at most m{n — s) + s < m min(n, to — 1) + rt. 

□ 

The parts of the algorithms that concern monomial manipulations are harder 
to analyze. As stated in the introduction, the number of integer comparisons 
in the monomial manipulation-part is reported to be proportional to n^TO^, 
assuming a standard order. We agree that this is an upper bound for the 
algorithm using steps C2' and C5'. Since it is not explained in [Ij how to check 
if t is a multiple of an element in L or of in(G), the actual behavior of this 
algorithm might be worse than the algorithm using step C2' and C5', although 
the list L during the former algorithm contains less elements than the list L 
during the latter. 

To check if i is a multiple of an element in in(G) using step C2' and C5' is 
simple and is due to the following nice observation given in [4 . 

Lemma 3.3. Using steps C2' and C5', to check if t is a multiple o/in(G) can 
be replaced by checking if |supp(t)| > cp(i), where cp{t) denote the number of 
copies of t in L. 

Proof. Let t be the first element in L. There exists exactly supp(t) submono- 
mials of t of degree \t\ — 1. All copies of t comes from a submonomial of of t 
of degree \t\ — 1 in B, and each submonomial gives rise to exactly one copy. It 
follows that if there exists a submonomial of t not in B then | supp(t)| > cp(t) 
and vice versa. But the first statement is equivalent with t being a multiple of 
an element in in(G). □ 
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We remark that this simphfication of the check does not apply when using 
steps C2 and C5. 

3.2 Optimizing the BM-algorithm 

In order to improve the time complexity during the monomial manipulations 
part, we will use a projection technique together with Theorem 12.151 The idea 
is to identify a set Ess of the variables with the property that supp(i3) C Ess 
and I Ess I < min(m — l,n). Once the set Ess is identified, we only consider 
monomials in the monoid generated by Ess, hence the associated order vectors 
for the monomials is of length bounded by min(TO — l,n). 

The projection technique is covered in the following lemma and is in some 
sense inspired by O [7] . 

Lemma 3.4. Let ^ he a monomial order and let pi, . . . be distinct points 
in k". Then we can determine a subset Ess of {xi, . . . ,Xn} with the following 
properties 

• supp(i?) C Ess 

• I Ess I < min(rn — l,n) 

• dimij(Ess(P)) = | Ess | ( Ess{P) has full rang.) 

e I{P) where Xk >- Xi^ if Ckj ^ 0. 

in 0{nnn?) arithmetic operations plus time complexity Varord^(n). 

Proof. We first determine ii, . . . such that Xi^ )~ ■ ■ ■ >- Xi^ which has time 
complexity Varord^(n). Let 



and let Ess = En. It is clear that if Xk ^ Ess, then Xk G in(/), since Xk can be 
written as a linear combination of smaller elements, hence supp(i?) C Ess. By 
a dimension argument, we have | Ess | < m — 1 and since Ess is a subset of the 
variables, clearly | Ess | < n. By construction dimi;(Ess(P)) — \ Ess |. To deter- 
mine if Xi^_-^^{P) ^ span]]5{l(P), £'j_i(P) is 0{m?) arithmetic operations by 
using a matrix representation. This is repeated n times, which gives the arith- 
metic complexity Olnm?). Finally when Xk ^ Ess we obtain a the expression 
Xk = 'Ylij c-kjXi. mod I{P) as a side effect of the matrix representation. □ 

Let TT be the natural projection from k" to k" with respect to Ess = 
{xij , . . . , Xi-}, that is, 7r((ai, . . . , a„)) = (a^j , . . . , a.;-). Let tt* be the corre- 
sponding monomorphism from T = h\yi-^ , • • ■ , Vi—] to S given by j/i . ^ Xi - . If f 
is any element in T, then by construction ■n*{f){p) — f{Ti{p)). This shows that 
if Ess(P) has full rang, so has (7r*)~^(Ess)(7r(P)). It follows that the points in 
7r(P) are distinct. This leads us to use the following isomorphism result. 



E. 




if j = 

} ifx,„_^.(P)^spanJl(P),i?,_i(P)} 



otherwise 
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Lemma 3.5. Let pi,...,p,„ be distinct points i k". Let I be the vanishing 
ideal with respect to these points. Let n be a projection from k" to k" such that 
7r(pi), . . . , 7r(pm) are distinct. Let T — k[yij , . . . yi—] and let J be the vanishing 
ideal with respect to 7r(pi), . . . , 7r(p„i). Then S/I and T/J are isomorphic as 
algebras. 

Proof. Let tt* be the map defined above. For / £ T we have 7r*(/)(p) = /(7r(p)). 
Notice that / € J is equivalent to f{n{qi)) = 0,Vi, which is equivalent to 
7r*(/)(pi) = 0,Vi, which is equivalent to 7r*(/) G /. This allows us to extend tt* 
to a monomorphism from T/J to S/I. Since 7r(pi), . . . , 7r(pm) are distinct, we 
have dimij(r/J) = dimi;(S'//) and thus, the extension of vr* is an isomorphism 
of algebras. □ 

Lemma 3.6. Let ■< be a monomial order on S . Let Ess be a subset of the 
variables such that supp(_B) C Ess, | Ess | < min(m — 1, n) and dimi;(Ess(P)) — 
I Ess I, where B is the monomials outside in(/(P)) with respect to -<. Let tt 
be the projection defined by 7r((ai, . . . , a„)) = (a^^ , . . . , ai_) and let tt* be the 
corresponding monomorphism from T ~ k[?/ij , . . . , to S. Let -<' be the 
monomial order defined by y" -<' if 7r*(y") -< TT*{y^). Let B' = {b[, . . . , 
be the set of monomials outside in(/(7r(P)) with respect to -<' . Then Tr*[B') = B. 

Proof. Suppose that 7r*(6^) can be written as a linear combination of elements 
in B; TT*{b[) = Y^jCjbj with ^T*{h'^) >- h, if ^ 0. Since supp(B) C Ess, 
we get b[ = Y^jCj{-K*)~^{bj) with 6- ^' {■K*)~^{bj), which is a contradiction. 
Hence tt*{B') C B from which it follows that tt*{B') — B hy Lemma [33] and a 
dimension argument. □ 

Lemma 3.7. In the context of Lemma \3.4\ and Lemma \3.6\. suppose that Xk — 
y^^. CkjXi- e I{P) for all Xk outside Ess. Let G' be a reduced Grobner basis for 
/(7r(P)).' Then 

G = 7r*(G') U {xk - ^Cfe,^*(Nf((^*)-i(x,J,G"))} 

j 

is a reduced Grobner basis for I{P). 

Proof. Since tt*{B') is the complement of in(/(P)) by Lemma [3.61 it follows 
that in(/(P)) is minimally generated by Ess'^ U7r*(in(/(7r(P)))). Clearly 7r*(G') 
is contained in G. Thus, it is enough to prove that 

Xk-J2ckj7T*{m((7r*y\x,^),G'))eI{P) 
i 

and that Xk is larger than any monomial occurring in the right hand sum. Since 
Xk - I]j Ck]Xi- £ I{P), we have that 

j 

is equivalent to 

j J 
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Using the monomorphism (tt*) ^ we see that this is equivalent to 

^c,,((7r*)-i(:.g-Nf((^*)-i(xg,G')e/(7r(P)). 

j 

Since each term {7T*)^^{xi. ) — Nf ((tt*)^^ (xi^. ), G") is in G' , we are done with the 
first part. The second part foUows since Xk >- Xi^ if Ckj ^ and each T:*{xi-) 
is written as a linear combination of elements less than tt* [xi . ) with respect to 

' □ 

We now give an example of our method. 

Example 3. Consider the points 

Pi = (1,1,0,1,0),P2 = (2, 2, 1,1,1), P3 = (2, 0,1, 1,-1), P4= (5,3,4,1,2) 

in Q'^, with respect to pure lex and xi )~ X2 )~ X'i )^ x/^ x^. Using Lemma \3.4\ 
we get Ess — {x^, x^} and 

X4 = 1,X2 ^ X5 + 1, Xi = X3 + 1, 

everything mod I{P). Thus, let 7r(ai, . . . , a5) = (03,05) and let T — k[2/i,y2] 
and let tt* he defined by yi t-^ X3 and 2/2 ^ ^5- We have that -<' is lex with 
Vi y' 2/2 and 

7r(F) = {(0,0),(l,l),(l,-l),(4,2)}. 

A call to the BM-algorithm with Tr{P) and -<' yields B' — {1,1/2,2/2,2/2} as the 
set of monomials outside in(/(7r(P))) and 

yt + 22/2 ~yl~ 22/1,2/1 - 2/i 

as a Grobner basis G' for I{'k{P)). Thus, a Grobner basis G for I{P) is 

{xl + 2x5 - 2:5 - 2x1, X3 - 

X4-7r*(Nf((^*)-i(l),G')), 
X2 - 7r*mi7T*)-\x5), G')) - ^*(Nf((7r*)-i(l), G')), 
x^ - 7r*(M{{n*)-\x3), G')) - 7r*(Nf ((^*)-i(l), G'))} 
and the complement o/in(/(P)) is 

B = {(7r*)-i(l), (7r*)-i(2/2), {^*)-\yl), {i^*)-\yl)} = {l,X5,xl,xl}. 
We have 

7r*(Nf ((7r*)-i(l), G')) = 7r*(Nf (1, G')) = 7r*(l) = 1, 
7r*(Nf((7r*)-i(x5),G') = ^*(Nf(2/2,G')) = 7r*(2/2) = X5 

and 

7r*Nf((7r*)-i(x3),G')) = 7r*(Nf(2/i,G')) = 7r*(2/|) = 

Thus 

G = {xl + 2x5 - 2^5 - '^xl,X3 ~ xl,X4 - 1,X2 - X5 - l,xi -xl- 1}. 
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Notice that although Ess ~ {a;3,a;5} were hnearly independent with respect 
to F, it did not foUow that Ess C B. 

To analyze the complexity of the method above, we first determine the cost 
of the monomial manipulations. 

Proposition 3.8. Let ~< be a standard order. An upper bound for the time 
complexity of the monomial manipulation part of the BM- algorithm using the 
projection technique is 

0(min(m, n)m^ log(m)) 

Proof. Suppose that m < n. We can now use the projection technique described 
above so that n < m — 1. Since the projection technique only affects the 
arithmetic complexity, we do not need to consider the cost for it in this analysis. 
By Lemma |2.6[ to determine ~<ess is 0(1). Everytime we insert an element into 
B, we insert at most n elements into L. Thus, the number of elements in L 
is bounded by nm. Thus, the complexity of the monomial manipulation part 
is dominated by merging a sorted list of n monomials with a sorted list of at 
most nm monomials, repeated m times. To compute ov{xiX°') given ov(a;") is 
0(log(m)) by Lemma Thus, each time an clement is inserted into B, it is an 
0(n log(m))-operation to create the list of monomials which we will merge with 
L. Since we create at most m such lists, the total time needed for creation is 
0{nm log(m)). Using Theorem l2.15l we see that each merge has time complexity 
0(jT^ max(log(m/'n), 1) + nmlog(n + m)). 

If TO > n, then all arguments hold if we replace n by n. Thus, in general, 
each merge has time complexity 

0(min(m, n)^ max(log(TO/ min(m, n)), 1) + min(TO, n)TOlog(min(m, n) + to)) 

which equals 

0(min(TO, n)m log(TO)) 

by a straightforward calculation. Since there are exactly m — 1 merges, we get 
the complexity 

0(min(TO, n)m? log(TO) + min(m, n)m log(m)), 
where the last term comes from the creation process and is negligible. □ 

Proposition 3.9. Let -< be an order defined by an integer matrix A. Let c — 
max(|aij|). We give two upper bounds for the time complexity of the monomial 
manipulation part of the BM-algorithm using the projection technique, based on 
two different methods. When m > n, the methods agree and an upper bound is 

0{n^m log(cTO) + nrn^ log(n + to) + ti^to log(cTO) + log(c) log(n)). 

When m < n, the first method has the bound 

0{m^ log(cTO) + log(c) log(7i)) 

to which one needs to add the cost for 0{nm^) arithmetic operations over Q. 
When m < n, the second method has the bound 

0{nm? log(cTO) + log(c) log(n)). 
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Proof. First of all we need to determine Varord^, an log(c) log(n)-operation 
by Lemma 12.101 Suppose that m < n. We can use the projection technique 
described above and we now have two choices. Either we use Lemma 12.111 to 
construct an n x n-matrix Aess using 0(nm?) arithmetic operations over Q, or 
we can use the n x n-submatrix of A, where we keep the columns that refers to 
the variables in Ess. 

In the first case, the cost for computing ov(xiTO) given ov(to) is 0(nlog(cm)) 
by Lemma [2.81 so the total time needed for the construction of the associated 
order vectors is 0(n^mlog(cm)). By Lemma and Theorem 12 . 151 we see that 
each merge has time complexity 

0(n^ log(cm) + nm log(n + m)) — 0{inn? log(cm)) 

as we merge a list of n elements with at most nm elements. Since we make 
TO — 1 merges, we deduce that the overall time complexity of the first method is 

0{m^ log(cm) + n'^m log(cTO) + log(c) log(n)) 

= 0(to'^ log(cTO) + log(c) log(n)) 

to which we need to add 0{nm?) arithmetic operations over Q. 

The second method differs from the first in that the vectors are n-tuples 
rather than n-tuples. Thus, computing ov{xi'm) given ov(to) is an nlog(cTO) op- 
eration, so the total time needed for the construction process is 0(nnTOlog(cTO)). 
Each merge requires time 0(nnlog(cTO) + nmlog(n-t-TO)) = 0(nm log(cTO)), so 
the overall time complexity of the second method becomes 

0{nrn? log(cTO) -|- nnTOlog(cTO) -|- n? log(c) log(n)) 

= 0{nrn? log(cTO) + n? log(c) log(n)). 

When m > n,we do not need to project and the two methods agree. The cost 
for the construction of the associated order vectors becomes 0(n^TOlog(cTO)), 
each merge is ri^ log(cTO) -I- nm log(n + m) and thus an upper bound is 

0(n^mlog(cTO) -I- nm^ log(n + to) + n^TOlog(cTO) + n^ log(c) log(n)). 

□ 

We are ready to state the main theorem. 

Theorem 3.10. An upper bound for the arithmetic complexity of the BM- 
algorithm using steps CI, C2', C3, C4 and C5' and the projection technique 
based on Lemma \3.4\ is 

0{nm^ + min(TO, n)m^). 
To this we need to add the time complexity 

0(min(TO, n)m^ log(m)) 

when ~< is standard. 

When -< is given by a matrix there are two methods to use. When m > n, 
the methods agree and an upper bound is 

0{n^m log(cTO) + nm^ log(n + to) -|- n^m log(cTO) -|- log(c) log(n)). 
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When m < n, the first method has the bound 

0{m^ log(cm) + log(c) log(7i)) 

to which one needs to add the cost for 0{nm?) arithmetic operations over Q. 
When m <n, the second method has the hound 

0{nm^ log(cTO) + log(c) log(n)). 

Proof. The complexity result for the monomial manipulation part is given in 
Proposition 13.21 If m < n and we use the projection technique described in 
Lemma 13.41 we will get the set B' of monomials outside in(/(7r(P))) and a 
Grobner basis G' for in(/(7r(P))). By Lemma [3.61 B — tt*{B') and by Lemma 

G = TT*{G') U {xk -J2ckj7T*{m{{n*)-\x,^), G'))}. 
j 

The computation of i^*{B) is immediate, since the only computation needed 
is changing indices. The arithmetic complexity bound follows if we can show 
that each Xk — '}2,j c^j tt* (Nf ((7r*)~^(a;i^. ), G')) is computable within 0{m?) arith- 
metic operations. To get a short proof, we will not use the information Xk — 
X]j '^kjXi^ G /. Instead we compute the evaluation vector {xk{pi), ■ ■ ■ ,Xk{pm)) 
= (pi/c, ■ • • ,Pmk) and write it as a linear combination of the elements in B, an 
operation which requires O(m^) arithmetic operations. Since i? is a basis, the 
linear combination will then equal J^j Cfej7r*(Nf ((7r*)~^(a;i^ ), G')). 

The complexity for the monomial manipulation part follows from Proposition 
13.81 and Proposition 13.91 

□ 

The following corollary states that our version of the BM-algorithm is pre- 
farable to the EssGB-algorithm [7]. 

Corollary 3.11. When m < n and the order is standard, the BM-algorithm 
using steps CI, C2', C3, C4 and C5' and the projection technique based on 
Lemma \3.4\ has arithmetic complexity 

To this we need to add the time complexity 

0{m^ log(TO)). 

3.3 Applications to the FGLM-algorithm and for ideals 
defined by functionals 

As was noticed in [TO], both the FGLM- and the BM-algorithm are instances 
of definitions of an ideal defined by means of a finite set of functionals Li : 
k[xi, . . . , Xn] k, such that / is in the kernel of \1/ : k[xi, . . . , Xn] k™, ^(/) — 
Li{f), . . . , L„i{f). For the BM-setting, the functionals are defined by Li{f) — 
f{pi) and in the FGLM-setting, the functionals are defined by Nf(/, Gi) = 
^ Li{f)ei. In [ro|, a list of different problems that can be seen as instances of 
an ideal defined by functionals is given. 
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If we use the steps C1,C2',F3,C4,C5' of the BM-algorithm where F3 is de- 
fined below, we obtain Algorithm 1 in jTU]. 

F3 Compute '^{t) — (61, . . . ,bm) and reduce it against the rows of C to 
obtain 

. . .,Vm) = (&i, ... ,6m) - ^aj(cji, . . . ,Ci„) flj e k. 

i 

It is reported in TO] Theorem 5.1, that Algorithm 1 in [10] needs 0{nm^ + 
fnrn^) arithmetic operations, where / denotes the cost of evaluating a func- 
tional. However, as for the BM-algorithm, one can replace the term nra^ by 

-|- nm? . Since we have shown (Proposition 13. 2p that the number of calls to 
C3 equals \G\ -f m = n -f min(m, n)m m, one can replace the term fnrn? by 
fnm + / min(TO, njm? + fm. Thus, Algorithm 1 in [lOJ uses 0(min(m, n)rrv^ + 
nm? -|- fnm -f / min(m, n)m?) arithmetic operations. In the BM- or the FGLM- 
situation, it is shown in [lO that / = 1, by a recursive argument. 

To the arithmetic complexity one needs to add the cost for the monomial 
manipulations, which is the same as for the BM-algorithm, since it is clear that 
we can use the projection technique described in Lemma 13.41 if we replace 

^spanJl(P),i?,_i(P)} 

by 

^{x,^-,) i spanj*(l), . . . , ^{x,J\. 

It follows that Theorem 13.101 can be lifted to the general setting of ideals 
defined by functionals. We state 

Theorem 3.12. An upper bound J or the arithmetic complexity of Algorithm 1 
in JlOjl using the projection technique is 

0(min(m, n)nC' + mri^ + fnm + f min(TO, n)m'^). 

To this we need to add the time complexity 

0(min(m, n)m^ log(m)) 

when ~< is standard. 

When -< is given by a matrix there are two methods to use. When m > n, 
the methods agree and an upper bound is 

0{n^m\og{cm) + nm^ log(n + rn) + n^mlog(cTO) -|- log(c) log(n)). 

When m < n, the first method has the bound 

0{m^ log(cm) -I- n^ log(c) log(n)) 

to which one needs to add the cost for 0{nm^) arithmetic operations over Q. 
When m < n, the second method has the bound 

0{nm^ log(cTO) -I- log(c) log(n)). 
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4 Discussion and future work 



The projection idea can be used for ideals defined by projective points and also 
in the non-commutative versions of the FGLM-algorithm, given in [2J . It should 
be remarked that the monomial manipulations with respect to non-commuting 
variables is computationally harder than those for commuting variables. 

Future work involves the question: Why do we need a Grobner basis for 
ideals defined by vanishing points? In the biological applications, where one is 
primary interested in normal form computations, it seems enough to compute 
a set B such that [B] is a k-basis for S/I. 
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